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ABSTRACT 

A two-dimensional electrodynamical model is used to study particle acceler- 
ation in the outer magnetosphere of a pulsar. The charge depletion from the 
Goldreich- Julian charge density causes a large electric field along the mag- 
netic field lines. The charge particles are accelerated by the electric field and 
emit 7-rays via the curvature process. Some of the emitted 7-rays may collide 
with X-ray photons to make new pairs, which are accelerated again on the 
different field lines in the gap and proceed similar processes. We simulate the 
pair creation cascade in the meridional plane using the pair creation mean-free 
path, in which the X-ray photon number density is proportional to inverse 
square of radial distance. With the space charge density determined by the 
pair creation simulation, we solve the electric structure of the outer gap in 
the meridional plane and calculate the curvature spectrum. 

We investigate in detail relation between the spectrum and total current, 
which is carried by the particles produced in the gap and/or injected at the 
boundaries of the gap. We demonstrate that the hardness of the spectrum is 
strongly controlled by the current carriers. Especially, the spectrum sharply 
softens if we assume a larger particle injection at the outer boundary of the 
outer gap. This is because the mean-free path of the pair creation of the 
inwardly propagating 7-ray photons is much shorter than the light radius so 
that the many pairs are produced in the gap to quench the outer gap. 

Because the two-dimensional model can link both gap width along the 
magnetic field line and trans-field thickness with the spectral cut-off energy 
and flux, we can diagnose both the current through the gap and inclination 
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angle between the rotational and magnetie axes. We apply the theory to the 
Vela pulsar. By eomparing the results with the EGRET data, we rule out 
any cases that have a large particle injection at the outer boundary. We also 
suggest the inclination angle of ainc ^65°. The present model predicts the 
outer gap starting from near the conventional null charge surface for the Vela 
pulsar. 

Key words: gamma-rays: theory -pulsars: individual (Vela pulsar). 

1 INTRODUCTION 

Rapidly rotating young pulsars such as the Crab and Vela pulsars emit high-energy 7-rays 
and are among the brightest sources in the 7-ray sky. In the pulsar magnetosphere, electrons 
and/or positrons are accelerated above 10 TeV. Because the pulsar is an electric dynamo with 
an available potential drop, ~ 1.3 x 10^^P~^/^P^:(5 Volt, where P is the rotational period 
and P_i5 is the time derivative of the rotational period in units of 10~^^ s/s, theoretically, it 
is suggested that the pulsars utilize some fraction of $a to accelerate the particles. EGRET 
instrument observed 7-ray emissions from 7 pulsars (Thompson et al. 1999). The observed 
modulation of light curves and spectral cut-off energy test the theoretical models; e.g., the 
polar cap accelerator model (Sturrock 1971; Ruderman & Sutherland 1975; Scharlemann, 
Arons & Fawley 1978; Daugherty & Harding 1996), the outer gap accelerator model (Cheng, 
Ho & Ruderman 1986a,b, hereafter CHR86; Romani 1996; Zhang & Cheng 1997; Hirotani & 
Shibata 1999a,b), and the inner annual gap accelerator model (Qiao et al. 2004). All models 
predict the particle acceleration in a charge depletion region called the gap, in which the 
charge density differs from the Goldreich- Julian charge density, with an electric field parallel 
to the local magnetic field hues. 

The outer gap model has been successful in explaining the observed main features of 
the 7-ray light curves, such as the two peaks in a period and the presence of the emissions 
between the peaks. CHR86 first investigated electric structure of the outer gap. By solving 
the Poisson equation for the vacuum gap, CHR86 showed that the inner boundary of the 
vacuum gap locates close to the null charge surface of the Goldreich- Julian charge density, 
and concluded that the outer gap is lying above the last open lines and is extending between 
the null charge surface and the light cylinder, the axial distance of which is zuic — c/ft, 
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where Q is the rotational angular frequency and c the speed of light. The last open lines are 
defined by the magnetic field lines tangent to the light cylinder. With CHR86's vacuum gap 
geometry, the two peaks in a period have been interpreted as an effect of aberration and 
time delay of the emitted photons (Romani & Yadigaroglu 1995, hereafter RY95; Cheng, 
Ruderman & Zhang 2000) . RY95 also reproduced phase separation of the pulses in different 
photon energies from radio to 7. 

Hirotani & Shibata (1999a,b, hereafter HS99) have solved electrodynamics in the non- 
vacuum outer gap (see also Hirotani & Shibata 2001a,b). Although HS99 worked in only 
one dimension along the last open line, they first connected consistently the field-aligned 
electric field in the outer gap with the curvature radiation and pair creation processes. The 
model reproduced the EGRET phase-averaged spectra of the Vela hke pulsars. 

With these successful works, the outer gap has been considered as a promising accel- 
eration model accounting for the high energy emissions of the pulsars. However, there are 
some problems in the previous works. For example, CHR86 discussed the electrodynamical 
structure of the vacuum gap, although they assumed an exponential grow of the number 
density in the trans-field direction. Such a large number of particles must partially screen the 
field-ahgned electric field on the magnetic field lines, on which they fiow. Another problem 
is the one-dimensionahty in the HS99 arguments. The one-dimensional model neglected im- 
portant trans-field effects caused by the curvature of the field lines. The 7-rays are emitted 
tangent to the local magnetic field lines, converts into pairs on different field lines, and as a 
result, makes a trans-field distribution of the current. It is obvious that the electrodynamics 
such as the strength of the field-ahgned electric field is affected by the trans-field structures. 
Furthermore, the one-dimensional electrodynamical model has used the trans-field thickness 
as a parameter to be determined by the observed fiux (Takata, Shibata & Hirotani 2004a, 
hereafter Paperl). 

Recently, the emission region extending to both the stellar surface and the light cylinder 
has been proposed to explain the presence of the observed outer-wing emissions and of the 
off-pulse emissions of the Crab pulsar (Dyks & Rudak 2003; Dyks, Harding & Rudak 2004). 
Mushmov & Harding (2004) proposed the slot gap accelerator model (Arons 1983), which 
bases on the polar cap model, as an origin of the emission region by Dyks & Rudak (2003). 
On the other hand, it has been suggested that the outer gap may have an elongated emission 
regions (Takata, Shibata & Hirotani 2004b, hereafter Paper2). Thus, the origin of the 7-ray 
emissions in the pulsar magnetosphere has not been conclusive up to now. 



4 J. Takata, S. Shihata, K. Hirotani & H.-K. Chang 

In Paper2, we have solved the electrodynamics in the outer gap with a two-dimensional 
model in the meridional plane. We extended HS99's one-dimensional model into two-dimensional 
one with the trans-field structure. This two-dimensional electrodynamical model links the 
one-dimensional electrodynamical outer gap model by HS99 and the three-dimensional ge- 
ometrical outer gap model by RY95. 

In Paper2, we showed that the cusp of the inner boundary is located at the position 
where the local charge density caused by the current is equal to the local Goldreich- Julian 
value. This implies that the inner boundary of the outer gap shifts toward stellar surface 
from the conventional null charge surface as the current increases. In Paper2, the elongated 
gap was predicted by the cases of the particle injection at the outer boundary of the gap. 
Thus, the outer gap indicates the acceleration region predicted by Dyks & Rudak (2003). 
However, because the spectrum was not investigated in Paper2, it was not certain that what 
kind of gap structure reproduces the observed 7-ray spectrum. 

The purpose of this paper is to calculate the 7-ray spectrum using the two-dimensional 
electrodynamical model proposed in Paper2 and to examine dependence of the spectrum on 
the current through the gap and the inclination angle between the rotational and magnetic 
axes. We can estimate fiux as seen from the Earth more correctly than the one-dimensional 
model, because we solve the gap structure in both longitudinal and trans-field directions. 
We discuss the particle acceleration in the pulsar magnetosphere, and constrain the model 
parameters in terms of the spectral cut-off energy and the flux. To compare with the obser- 
vations, we relax the previous restrictions of the model such as 

• the inclination angle is equal to zero (i.e. aligned-rotator), 

• X-ray fleld of the pair creation is homogeneous in the magnetosphere, 

• all 7-rays produced by the curvature process have critical energy. 

In this paper, we calculate the pair creation process with inhomogeneous and anisotropic 
X-ray fleld and the curvature spectrum of the individual particles (§ ^2.4j) . 

We demonstrate that the calculated spectrum softens if one assumes a larger current 
through the gap. Especially, the calculated spectrum sharply softens with the increasing the 
number of particles injected at the outer boundary of the outer gap. On the other hand, we 
also demonstrate that an increase in the number of particles injected at the inner boundary 
or of the inclination angle moderately changes the calculated spectral hardness. Using this 
dependence of the spectral hardness, therefore, we can investigate the current in the outer 
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gap of the observed 7-ray pulsars by comparing the model spectrum with the observed one. 
We apply the theory to the Vela pulsar, and show the gap geometry, the current through 
the gap and the inclination angle (= amc) in the Vela magnetosphere. We show that the 
EGRET spectrum predicts a large inclination angle ^ 65° and only a few particle 
injections into the outer gap at the outer boundary. 

In ^ we present the basic equations describing the stationary outer gap structure. In 
^ ^2.41 we develop the previous pair creation cascade model in Paper2 using the Monte Carlo 
method. In §21 we compare the results with the observation of the Vela pulsar. In this paper, 
the polarity is f2 ■ m ^ 0, where m is the magnetic moment of the star. Since we ignore the 
effects of ions, the results do not depend on the polarity. 

2 MODEL & BASIC EQUATIONS 

For solving a stationary structure of the outer gap, we treat the ideal conditions: 

Caul the stationary condition {dt + ^Id^ = 0, where t and are the time and azimuth, 
respectively) is satisfied, that is, there is no time variation of any quantities as seen in a 
co-rotating system, 

Con2 the star crust is a rigid rotating perfect conductor, 

Con'd) the plasmas fill the pulsar magnetosphere, and the condition E • B = holds all 
the way along the field lines connecting the star and the outer gap accelerator. 

Strictly speaking, the condition Con3 cannot arise in the first place. We assume that a 
mechanism of the electric field screening works in the region between the stellar surface and 
the inner boundary of the outer gap. In S ^4.3[ we discuss the dynamics outside the gap. 

By neglecting the magnetic component by rotational effects and by the current, we 
assume the static dipole magnetic field. This static dipole structure will be a good approxi- 
mation as long as the electrodynamics in the outer gap is discussed. In fact, the sweepback 
effect is important in the geometrical studies to produce the double peaks in the light curves 
(Cheng, Ruderman & Zhang 2000). Electrodynamical speaking, the magnetic field affects 
on the electrodynamics thorough the Goldreich- Julian charge density, that is, the magnetic 
component projected to the rotational axis. As compared with the static dipole field, the 
null charge surface of the Goldreich- Julian charge density on the last-open lines, and there- 
fore the inner boundary of the outer gap for the rotating dipole field is located closer to the 
stellar surface about 10% of the light radius. With this small difference of the positions of 
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the inner boundaries between two magnetic fields, the quahtative arguments on the electro- 
dynamics in the outer gap will not be changed very much as long as the outer boundary of 
the gap is not located near the light cylinder. In the later sections, we will predict that the 
out boundary that is favored for the Vela spectrum does not locate near the light cylinder 
(FigE)), while the outer boundary on the light cylinder has been assumed in the previous 
light curve models (Romani 1996; Cheng, Ruderman & Zhang 2000). It is also worth noting 
that even if the outer boundary is located inside of the light cylinder, the two peaks in 
the light curves are expected because the outward moving particles accelerated in the gap 
continue to emit the high energy photons outside the gap via the curvature process (S §4.2j) . 

2.1 Poisson equation 

With the condition Conl, the electric field can be written as (Shibata 1995, Mestel 1999) 

E(r, t) = -(fie, X r) X B/c - V<l>„,o(r, t), (1) 

where $„co is the non-corotational potential and e, is the unit vector along the rotation 
axis. The E'n -acceleration is indicated by the non-corotational part in equation (^, E\\ = 
— B ■ V^nco/B 7^ 0. Substituting equation ((H) in V ■ E = Anp gives the Poisson equation, 

A$„,,o(r) = -47r[p(r) - pGj{r)], (2) 

where p(r) is the space charge density, and Pgj{^) is the Goldreich- Julian charge density. 



PGJ 



B - X (V X B) 



(3) 

To simplify the geometry, we assume that the gap dimension in the azimuthal direction 
is much larger than that in the meridional plane. Neglecting variation in the azimuthal 
direction, we rewrite equation Q as 

A,,e$„,o(r) = -47r[p(r) - PGj(r)], (4) 

where /\r,e represents (r, 6')-parts of the Laplacian. 

At the moment, let's rewrite down the Poisson equation (j3]) as 



dsl +^--Mp-Pc.), (5) 

where S|| and s± are the distance along the magnetic field and trans-field directions, respec- 
tively. Near the boundaries of the gap or for a thick outer gap, in which typical gap width 
along the filed line is much larger than the typical trans-field thickness, the left hand side of 
equation (jSj) is dominated by the first term. The accelerating electric field is determined by 
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^11 = -^^^ ~ / Att{p - pGj)ds\\. (6) 



ds\\ 

Apart from the boundaries for the thin outer gap, on the other hand, the second term 
dominates in the left hand side of equation (jS)). Hence, the accelerating electric field behaves 
such as (CHR86) 

E\\^2-Ks^{s^-D^)^{p- Pgj), (7) 

where we ignore the trans-field dependence of the charge density, is the typical gap 
trans-field thickness, and = and s_[_ = Dj^ represent the lower and upper boundaries 
of the gap, respectively. We find that the magnitude of the gradient of the effective charge 
density {= p — pcj) along the field lines determines the strength of the field-aligned electric 
field for the thin gap geometry, although the magnitude of the effective charge density causes 
the field-aligned electric field for the thick gap geometry. In a real pulsar magnetosphere, 
the thin gap or the thick gap will be controlled by the pair-creation mean free path. 

To solve the Poisson equation numerically, we introduce an orthogonal curvilinear 
coordinate system (x, () following the dipole field on the meridional plane as 

X = — '■ ^ — ^ = constant along a dipole field line, (8) 

and 

C = = constant along a curved line (9) 

^Ic 

perpendicular to the field lines, 

where r and 9 is the radial distance and the colatitude angle with respect to the rotational 
axis, respectively. The Laplacian becomes 

_3cos2 0'+l92 4 d 3cos2 0' + l92 3(3 cos^ 0' + 1) 9 

sin^ 9' xsin^O'd^^ 4cos3^' ^ 4Ccos3 9' ^C' ^ ^ 

where 9' = 9 — ainc- Non-corotational electric field (= E^co) is given by 

d^nco _ V3 COS^ 9' + 1 d^nco 



d^nco _ V3 COS^ 9' + 1 d<!?nco 

^" ds^ ~ zui,sm'9' dx ' ^ ' 

where ds^ = (2 cos^/^ 9'/y/3 cos^ 9' + l)dC and dsi_ = (sin^ 9'/y/3 cos^ 9' + l)dx are the line 
elements along the field line and the perpendicular curved line, respectively. 
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2.2 Particle continuity equations and pair creation process 

For the young pulsar such as the Vela pulsar, the non-corotational drift motion of the charged 
particles due to the non-corotational electric field E„co is negligibly small in comparison with 
the corotational motion (Paper2). Also ignoring the gyro- motion, we denote the velocity v 

as 

■y ^v\\e\\-\-wQ.e^, (13) 

where f || represents the longitudinal velocity. Because the accelerated particles move at the 
speed of light, the longitudinal velocity in the meridional plane is to be v\\ 
with which the particles migrate along the line of x = constant. 
With Conl, the continuity equation yields 

B.v(^)=±SW. (14) 

where S{r) is the source term, and and N_ denote the number density of outwardly and 
inwardly moving particles (i.e. the positrons and the electrons), respectively. 

To calculate the source term S{r), we adopt the pair creation process between background 
X-ray photons and the 7-ray photons in the gap. The mean free path of a 7-ray photon 
with energy is written as 

with 

where dEx ■ dNx/dEx is the X-ray number density between energies Ex and Ex + dEx, 
Ox-y is the collision angle between an X-ray photon and a 7-ray photon, Efh — 2(meC^)^/ (1 — 
cos9x'y)E^ is the threshold X-ray energy for the pair creation, and ap is the pair creation 
cross-section, which is given by 



f°° dNx 

r]p{r, E^) = (1 - cos^x7)c / dEx—=^{r, Ex)crp{E^, Ex), 



3 

ap{E^,Ex) = Y^^t(1 
where 



V 



2\ 



(3_^;4)inl±^_2^;(2-^;2) 
1 — V 



(16) 



(mgC^)^ 



v{E^, Ex) = . - ^ ^ 

\ 1 - cos Ox-f EjEx 

and (Jr is the Thomson cross section. In the latter sections, we apply the theory to the 

Vela pulsar, for which the observed X-ray spectrum is dominated by thermal emission. This 

component has been interpreted as the emissions from the stellar surface. At the radial 



A two-dimensional electrodynamical study on the gamma-ray pulsar 9 
distance r from the centre of the star, the thermal photon number density between energy 
Ex and Ex + dEx is given by 

dNx _ 1 /2ny /Reff\' El ^^^^ 



dEx ATxKchJ V r / e:>cp{Ex/kTs) - V 
where Reff is the effective radius of the emitting region, and kTg refers to the surface 
temperature. For the values of i?e// and Tg, the observed ones are used. 

With the soft photons from the stellar surface, the coUision angle Ox-y of the 7-ray photon 
after traveling the distance s from the emission point (ro, ^0) is obtained from 

S + rCOS^em f^^. 

cos6'x7(r) = , (18) 

where Of^rn is the angle between the emission direction and the radial direction at the emission 



point, which are cos^^ = y 1 — {wqVl/ cYBr{ro, 9o)/B{ro, Oq) for the outwardly propagating 
7-rays, and cos^^ — —^1 — {wofl/cyBr{ro,9o)/B{ro,9o) for the inwardly propagating 7- 
rays. For ccj^c — 45°, for example, the 7-rays are emitted outward (or inward) from the 
conventional null charge point on the last open line with 9'^^ ~ 22.8° (or ~ 157.2°). 



2.3 Curvature radiation process 

We calculate the curvature radiation process of the accelerated particles. The power per unit 
energy emitted by the individual electrons (or positrons) is written as 

P,iR,,r,E,)^^^F{x), (19) 



"'C 



where x = E^/Ec, 
3 hcT^ 

and 

POO 

F{x) = X / K^/z{t)dt, (21) 

where Rc is the curvature radius of the magnetic field line, F is the Lorentz factor of the 
particles, ^5/3 is the modified Bessel function of the order 5/3, h is the Planck constant, 
and Ec gives the characteristic curvature photon energy. The Lorentz factor in the gap is 

given by assuming that the particle's motion immediately saturates in the balance between 
the electric and the radiation reaction forces, 

/od2 \1/4 

Vsat{Rc.E\\)={-^E\\ + l\ . (22) 
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As discussed in Paper 1, the curvature emissions of the outwardly moving particles from 
outside the gap also contribute to the total emissions if the outer boundary of the gap is not 
located close to the light cylinder. Outside the gap, the particles loose their energy via the 
curvature emissions so that 

We assume that the particles migrating on the magnetic field line labeled by x escape 
from the gap with a Lorentz factor given by the value at iy||(x)/c = t^, which describes 
the condition that the gap crossing time Vr||(x)/c of the particles is equal to the radiation 
damping time (= td) estimated by 



2 



t,^A-W-' ^ s. (24) 

VlOOrads-V VlOV \0.5wicJ ^ ' 

In the present paper, we take into account the curvature emissions up to = 0.9cj7;c. 

The saturation in the gap simplifies the problem significantly. The saturation mo- 
tion of the particles will be achieved if the the typical acceleration length scale lac ~ 
10'''(r/10'')(-E'll/3 ■ 10''V • m^^)^^ cm is shorter than the gap width (see Paperl). 

For some part of the gap, where the electric field is week, the saturation approximation will 
break down (Hirotani, Harding & Shibata 2003). In the two-dimensional model, specifically, 
around the lower and the upper boundaries, the condition W\\ > lac is not satisfied effectively 
for all particles. Nevertheless, we adopt the Lorentz factor Tsat to all particles in the outer 
gap for simplicity, because the saturation breaks down within a few percent of the gap thick- 
ness measured from the upper and the lower boundaries and because the contributions of 
the emissions from such region are less important for the observed spectrum above lOOMeV. 

The high-energy particles also emit high-energy photons via the synchrotron process and 
the inverse-Compton process. In ^ ^4.H we show that the both emission processes will be less 
important for the gap electrodynamics in comparison to the curvature process. 



2.4 Pair creation position 

To calculate the source function S{r) and 7-ray spectrum, we apply the simulation code 
developed in Paper2 (but we will revise the code for some effects as described later): Stepl) 
We calculate number (= Gq) of the emitted photons per unit time per unit volume at 
position Xo (see Fig^. The 7-rays are emitted in the direction of the particle motion v at 
Xq. Step2) We calculate photon number converted into pairs in a distance so(~ ^ic) by 
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.out 




.ight cylinder 



Figure 1. The pair creation cascade model in the meridional plane. The electrons or positrons are accelerated by the electric 
field along the magnetic field lines and emit the 7-rays. In the meridional plane, the 7-rays radiated by the curvature process 
are beamed in the direction of local magnetic field. The radiated 7-rays may convert into the pairs by the pair creation process. 
The new born pairs are also accelerated and emit the 7-rays. The Xh Xupi Cim ^^^d Cout present the lower, upper, inner, and 
outer boundaries of the gap, respectively. 



6G^ = Gq[1 — exp(— /q*" 1/lpds')]. Step3) We divide 6G^ photons into n flux elements, and 
pick up n points on the three-dimensional path length (sj < So,i = 1, ■ ■ -n) following Monte 
Carlo method f^ ^2.4.2|) . The pair creation point X, in the meridional plane is deduced from 



|Xj — Xo| = SjCos6'^, where 9^ = cos^^ yl — {Qzu/wicy. Step4) We calculate the number 
of pairs produced per unit time per unit volume as the source term 5'(Xj). 
In this paper, we take into account the following effects: 

• the energy spectrum of the curvature photons emitted by the individual particles, 

• the soft photon field proportional to the inverse square of the distance (that is, in- 
homogeneous) and the soft photons propagate radially from the stellar surface (that is, 
anisotropic), 

which were ignored in Paper2 for simplicity. In the following subsections, we detail how the 
above effects are built onto the cascade model. 

2.4-1 Curvature spectrum 

To give Gq in Step3) , we define the following 7-ray photon number densities in the dimen- 
sionless energy interval between Xi and Xj+i, 




n 



j::^'p,{x)dx 

{xi + Xi+i)/2 



v^e^r F{x)dx 
hRc {xi + Xi+i)/2 



(^ = 1,---, 



ni-1) 



(25) 
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where Xi = EljEc. We define xi = and select each value of Xj {i = 2, . . .rii) so as to 
satisfy j^^^^ F{x)dx = ai/rii, where Oi = F(x)dx = Stt/S^^^. Instead of x = oo, we 
define Xm+i = (3x„i — x„j_i)/2 and Gq^ = ai/xm+i- For each bin Gq, we carry out the 
processes Stepl)~Step4). In this paper, we adopt rii = 20 and n = 20 in Step3). 



2.4-2 Pair creation cascade model 

The radial dependence in the mean free path of equation p5|l can be written as Ip oc 
(1 — cos6x'y)^^r'^, where arises due to soft photon field diluting with the radial distance 
and 6'x7 of equation p8|) is the collision angle between an X-ray photon and a 7-ray photon. 
Because cos Ox-y ~ cos 9em is satisfied as long as the propagating distance s is smaller than 
the radial distance r ~ O.Stujc, the radial dependence of the collision angle for 7-rays emitted 
at the same position is less important. In this paper, therefore, we adopt 6x-y = Oem, which 
are collision angle at the emission point of the 7-rays considered. We take into account the 
fact that the value of depends on the emission point and the emitted direction (inwardly 
or outwardly) of the 7-rays. 

Using Monte Carlo method, we simulate the pair-creation position in the gap. The prob- 
ability that a 7-ray photon makes a pair between distance s and s + ds from the emission 
point is 

Ppds = ^"P(-/oy;^(^>^O rf,. (26) 

Using a function [= f{y)], we transform uniformly distributed random numbers, {yi,y2, • • • , Vn} = 
y, into random numbers, {si, S2, ■ ■ ■ , s„} = s, following the probability Pp. We obtain the 
function f{y) by solving the equation (Press et al., 1988) 



For the present case, /p oc r^, we can perform the integration in equation ()27|) analytically. 
As a result, we obtain the function f{y) as 



sin Oem^T^y , . f ^ 
-tq h atan 



Ci Vtan 6*, 



'em 



r COS 9em, (28) 



s = f{y) = ro sin 6*6™ tan 
where 

27v^ RLf El 
c, = (1 - -osex,)^^ j^^^^ '^^^ eMEx/m-l '''- 

The transformation represented by equation ()28|) gives a random number s of < s < cxo 
with a uniformly selected random number y of < ?/ < 1 . Because we are concerned with the 
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pair creation positions between < s < Sq ~ 'UJic, we define uniformly distributed random 
numbers as = (cj^ — l)y + 1, wliere 
f 1 



exp 



ro sin 9e 



^ (So + ro cos 9em\ ^ f I 

atan — atan 



r sin 9em J \tan6', 



A series of the transformations, y ^ ^ s = /(y^), produces the pair creation positions 
s following the probability Pp between < s < Sq ~ 

Following above Monte Carlo method, we choose n pair-creation positions for Gq photons 
emitted at the position considered. On each pair-creation position, the pairs produced per 
unit time and per unit volume is given by Ggfl — exp(— l/lpds')]/n. 



2.5 Boundary conditions 

We introduce the conventional dimensionless variables: 

(29) 



Up / 47re^ ^Bic 

c ' ^ \ nip 27rce' 



$ = $ E = -V<l> (30) 

J± = ^^, B = B/5,„ (31) 

/? = ^, (32) 
and 

P=nB^j2^c = ^,^^^'^-^^ PC. = -5.. (33) 
The Poisson equation ^ and the continuity equation (fT^ are rewritten as 

Mnco = -[w^U-J-) + B,l (34) 

PI I 

and 

B-Vj± = ±5(f), (35) 
respectively. 

In the meridional plane, the outer gap has four boundaries called as upper, lower, inner 
and outer boundaries (see FigUI). The upper and lower boundaries are defined by the mag- 
netic surfaces labeled by Xup and Xh respectively. The inner [= Cm(x)] and outer [= Coutix)] 
boundaries are functions of x- Fo^^ example, CoMt(x)=constant represents the outer boundary 
defined by a curve line perpendicular to the magnetic field lines. 
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We impose some conditions on the four boundaries. The inner and the outer boundaries 
are defined by the surfaces on which the field-ahgned electric field E\\ vanishes, 



With Con?) in §21 we have postulated that the condition E ■ B = holds between the stellar 
surface and the outer gap. The uniform value of $„co over the stellar surface propagates 
to the inner, upper and lower boundaries. By setting an arbitrary constant for ^nco on the 
stellar surface to be zero, we impose 



As discussed in § §4.3[ there is a possibility that the condition E • B = is not held between 
the stellar surface and the inner boundary. For such a case, it is difficult to determine the 
boundary condition on the inner boundary of the outer gap. However, because the potential 
drop between the stellar surface and the inner boundary is expected much smaller than the 
potential drop of whole outer gap as discussed in ^ ^4.31 the condition (jSZI), ^nco = 0, at the 
inner boundary is a good treatment. 

The continuity equation (jH^j) satisfies the conservation law of the longitudinal current 
density, that is, the total current density j+ + j_(= jtot) is constant along a field line. This 
total current jtot can be rewritten by sum of external and internal components. The external 
component is separated into two parts: 

(1) ji, the current carried by the positrons (e.g. originating in sparking on the stellar 
surface) coming into the gap through the inner boundary, 

(2) j2, the current carried by the electrons (e.g. originating in pulsar wind region) coming 
into the gap through the outer boundary. 

In the following, we call ji and j2 as the external positron and electron components, respec- 
tively. The internal component (= jg) is carried by the electrons and positrons produced in 
the gap. In terms of (jg, ji, ^2), the total current density becomes 



The current should be determined by some global conditions, because the outer gap, the 
polar cap and the pulsar wind interact each other through the current. By parameterizing 
external positron and electron components, we impose the conditions as follows 




(36) 




(37) 



(38) 




(39) 
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and 

J-{Cout)=j2{x)- (40) 

We have four model parameters, i.e. the inchnation angle (ai„c) and the three current 
components (ji, J2, Jg)- Because we have Dirichlet-type and Neumann-type conditions on the 
inner boundary, we solve the position of the inner boundary. We fix the lower boundary to 
the last open line. In the present model, we obtain unique positions of the inner and outer 
boundaries by giving the three current components and the position of the upper boundary. 
In the calculation, however, we give the position of the outer boundary instead of the internal 
current {jg) and solve the trans-field distribution of jg, because it is difficult to manage the 
distribution of jg, which is affected by the pair creation cascade, by hand. 

We consider the following particular solution in the present paper. In general, the inter- 
nal current jg increases as the outer boundary shifts toward the light cylinder if one fixes 
the position of the upper boundary. In the present electrodynamical model, the solutions 
satisfying the boundary conditions disappear if a large internal current jg fiows on a mag- 
netic field line. We regard that the gap having a marginal structure represents the state that 
actually appearers. This critical value of the internal current component is denoted by jcrit- 
If one gives the position of the upper boundary ixup) and the external current components 
(ji, J2)) the trans-field distribution of jg having the critical internal components jcrit and the 
positions of the inner {Cm) and outer {Cout) boundaries are determined. 

As we just mentioned, the solutions satisfying the boundary conditions disappear if the 
internal component jg becomes too large (> 0.1 ~ 0.2), because the produced pairs may 
make the super Goldreich- Julian charge density and may quench the gap. On the other hand, 
we can significantly increase the external component ji or j2 because the injected particles 
only shift the gap outward and inward. For example, if the external electron component 
increases, the effective charge neutral surface and the outer gap position shift inside relative 
to the conventional null charge surface of the Goldreich- Julian charge density (see § §3.111 . 

A solution is given by a self-consist connection between the field-aligned electric field 
determined by the distribution of the charge density (i.e. the pair creation cascade), the pair 
creation cascade controlled by the curvature radiation process (i.e. the Lorentz factor in the 
gap), and the Lorentz factor calculated from the electric field. 

The expected 7-ray fiux on the Earth is given by 

^(^7) = 7^^7(^7) (41) 
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where Ai7 is solid angle of the 7-ray beam and N^{E~^) is the total photon number emitted 
per unit energy per unit time, 

M^(E^)^j PJ^^W^iv, (42) 

where Vtot represents the total volume of the emission region, which includes both inside and 
outside the gap. In the present two-dimensional model, the solid angle of the 7-ray beam is 
estimated as AVL ~ 2txD ^jrout-, where and tout are, respectively, the solved trans-field 
thickness and the radial distance to the outer boundary of the outer gap. The total photon 
number N-y{E^) emitted in the gap depends on the azimuthal spread angle (= A0) of the 
outer gap, which is unsolved in the present model. Because there are no detail discussions 
for the pair creation cascade process with the three dimensional geometry up to now, it is 
difficult to estimate the azimuthal spread angle of the gap in the present two-dimensional 
model. Therefore, we use the azimuthal spread angle as a model parameter to be determined 
by the observed flux. The observed emission phase in the light curves (Kanback et al. 1994; 
Fierro et al. 1998) and the previous geometrical studies (Romani 1996; Cheng et al. 2000) 
for the light curves model have predicted ~ vr radian. 



3 RESULTS 

In this section, we present the results where we use the parameters of the Vela pulsar whose 
X-ray field at the outer gap is deduced from the observations. From ROSAT (Ogelman, 
Finley & Zimmerman 1993), RXTE (Harding et al. 2002) and Chandra (Pavlov et al. 2001) 
observations of the Vela pulsar, the spectrum in X-ray band is expressed by two components; 
power- low component and surface black body component. Because the former component 
is much fainter than the latter component, we neglect the former contribution to the X-ray 
field of the pair creation process. Ogelman et al.(1993) fitted the thermal component with 
T°° ~ 1.7 MK and Reff ~ 1.3((i/0.25kpc)^ km, where d is the distance to the Vela pulsar. 
In this paper, we adopt T°° = 1.7 MK and d = 0.25 kpc (Cha, Semback & Danksm 1999). 

In the following subsections, we examine the dependence of the 7-ray spectrum on the 
current components by fixing firstly the critical internal component jcrit (§ ^3.111 and secondly 
the external components, ji and j2 (§ ^3.2|1 . In § §3.3^ we compare the results with the observed 
spectrum of the Vela pulsar, and diagnose the outer gap geometry, the current and the 
inclination angle. 
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Figure 2. The trans-field structures of the internal component jg for case45l~case45 4. 

3.1 Gap structures for various external components 

In this subsection we fix tlie inclination angle at ainc = 45° and the critical internal compo- 
nent at jcrit ~ 0.12. To examine the effects of the external components (ji, 72); "we consider 
the following four cases: 

case45 1: (ji, J2) = (10"^ 10"^) for Xup>X> Xh 

case45 2: (ji,j2) = (0.05,0.05) for Xup>X> Xh 

case45 3: (ji,j2) = (0.1,0.0) for Xup>X> Xh 

case45 4: (ji,j2) = (0.0,0.1) for Xup>X> Xi- 

In case45 2~case45 4, the external current is 10% of the Goldreich- Julian value, which is 
carried by both components (for case45 2), the positron component (for case45 3), and the 
electron component (for case45 4). In case45 1, we consider a case of very small external 
component, that is, whole current is carried by the internal component, jtotix) ~ jgix)- 

FigOshows the trans-field distribution of the internal component jg for case45 l~case45 4. 
The abscissa and ordinate refer, respectively, the current density in units of the Goldreich- 
Julian value and the hight from the last open line (i.e. the lower boundary). The total current 
density on each magnetic field line is represented by jtotix) = jgix) + Ji + J2- So, for each 
case, the total current density on the magnetic field line where the internal component has 
the critical value {jcrit ~ 0.12) is jtot ~ 0.12 for case45 1 (ji = j2 = 0) and jtot ~ 0.22 for 
case45 2~case45 4 (ji + j2 = 0.1). 
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Figure 3. The geometries of the outer gap (filled region in each panel) for case45 l'~case45 4. The solid and dashed and 
dashed-dotted lines show the last open line, the curve line representing the outer boundary and the conventional null charge 
surface, respectively. On the dotted line in each panel, the condition, jcrit + 32 — jl = B^/B, where jcrit ~ 0.12, is satisfied . 
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Figure 4. Par-creation probability of a 7-ray with IGeV emitted from the null point on the last open line. The solid and 
dashed lines refer, respectively, the probability of the inwardly and outwardly propagating 7-rays in the Vela magnetosphere. 
The dashed line shows the value increased 10^ times of the probability. 



3.1.1 General features 

The filled regions in FigElshow the outer gap in the magnetosphere for case45 l~case45 4. 
The inner boundary in each panel satisfies both conditions E\\ = and $.„co = 0. The solid 
and dashed-dotted lines in each panel are the last open line and the conventional null charge 
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Figure 5. The pair creation position map in the gap of case45 2. The whiteness refers the particles flux (divided by f2_B/27re) 
produced per unit length [/cm] due to the pair creation at each position. 

surface, respectively. The outer boundaries of case45 1, 2, 3, and 4 are defined by the curved 
hues (the dashed hues) labeled by C^out ~ 0.9, ~0.65, ~0.95, and ~ 0.45, respectively. 
We list the general features that were already detailed in Paper2, in which aj„c=0°: 

• The external component affects the gap position in the pulsar magnetosphere. As ap- 
peared in Fig|2l the typical gap position of case45 3 or case45 4 is shifted outside or inside 
relative to the typical positions of case45 1 and case45 2. Specifically, if ji = j2 such as 
case45 1 and case45 2, the inner boundary of the lower parts of the gap, where jg ~ 0, is 
located close to the conventional null charge surface (pcj = 0). On the other hand, if ji ^ j2 
such as case45 3 or case45 4, the outer boundary shits outside or inside. This is because the 
external component shifts the effective charge neutral surface relative to the conventional 
null surface. 

• As the dotted lines, on which jcru + 32 — ji = B^/B with jcrit ~ 0.12 is satisfied, 
shows, the cusp of the inner boundary of each case is located at the position where the space 
charge density caused by the current is equal to the local Goldreich- Julian charge density. 
The inner boundary shifts toward the stellar surface as the critical and/or external electron 
components increases. 

• As FigEl shows, the present model predicts that more large current runs though the 
middle part of the gap rather than the upper part predicted in the CHR86 arguments. The 
internal component jg increases with hight in the lower half part of the gap, and decreases 
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in the upper half part. This decrease is because most 7-rays emitted in the lower part of 
the gap escape from the inner boundary or the outer boundary before reaching the upper 
part. Another reason is that the field-aligned electric field, the Lorentz factor and in turn 
the emitted photon energy decrease with increasing hight in the upper half part. 

3.1.2 Gap geometry & pair creation mean- free path 

In Paper2, the trans- field thickness of the outer gap was not sensitive to the external com- 
ponent (for example fig. 6 in Paper2), because we assumed a homogeneous X-ray field. On 
the other hand, we find in FigOlthat the trans-field thicknesses of case45 2 and case45 4 are 
much smaller than those of case45 1 and case45 3. This difference in the thickness is explained 
in terms of the mean-free path depending on the radial distance and the emission direction, 
which were ignored in Paper2, as follows. 

FigE] shows how the pair-creation probability of equation (j^UI) develops along the path 
of IGeV 7-rays emitted from the null point on the last open line. The solid and dashed lines 
refer the probabilities of a photon emitted inward and outward, respectively. The dashed 
line shows the value increased 10^ times of the probability. We find that the mean free path 
of the outwardly propagating 7-rays is much longer than that of the inwardly propagating 
ones, in which a photon out of about 10'^ photons converts into a pair after running about 
the light radius. FigEl shows the pair creation position map of case45 2 (ji = 72 = 0.05), 
where we do not plot the pairs created outside the gap. In Fig|3 the increase in the whiteness 
refers the pair flux produced per unit length at each position. As we expected, the inwardly 
propagating 7-rays make many pairs around the inner boundary and the middle part of the 
gap. 

The large difference in the mean-free path between the outwardly and inwardly propa- 
gating 7-rays is mainly attributed to difference in the collision angle Ox-y with the surface 
X-ray photon. For example, the collision angle of the 7-ray photon emitted from the conven- 
tional null point on the last open line of ainc = 45° becomes Ox-y ~ 22.8° for the outwardly 
propagating and becomes Ox-y ~ 157.2° for the inwardly propagating 7-rays. So, the pair 
creation processes occur with tail-on like collision (1 — cos6'x7 ~ 0) for the outwardly and 
head-on like collision (1 — cos^x7 ~ 2) for the inwardly propagating 7-rays. For the out- 
wardly propagating 7-rays with IGeV, the threshold energy of the X-ray photon of the 
pair creation is given by Eth = 2(meC^)^/(l — cos^x7)/lGeV ~ 6.7keV, which indicates the 
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Figure 6. The electric field structure for case45 l~case45 4. The solid, dotted, and dashed lines in each panel show the 
distribution of the field-aligned electric field on the three different magnetic field lines, which locate, respectively, 25, 50, and 
75% of the gap trans-field thickness measured from the last open field line. The abscissa represent the arc length with the origin 
at the null surface. 



deep Wien region of the Planck distribution for the Vela pulsar {kT^ ~ 0.15keV). For the 
inwardly propagating 7-rays, the threshold energy becomes Eth ~ 0.27keV. This difference 
in the threshold energy and the resultant difference in the X-ray photon number for the 
collision mainly produce the large difference in the mean free path. 

Because a large production of the pairs immediately quenches the gap, the mean-free 
path of the pair creation controls the gap size. As a results, case45 2 and case45 4, in which 
many 7-rays are emitted inward by the external electrons, have much smaller gap thickness 
than case45 3, in which many 7-rays are emitted outward by the external positrons. 



3.1.3 Electric structure 

Figini shows the distribution of the field-aligned electric field along the magnetic field lines 
for case45 l~case45 4. The solid, dotted and dashed lines represent, respectively, the field- 
aligned electric field on the magnetic field through the gap at 25, 50 and 75% of the trans-field 
thickness from the lower boundary. The abscissa refers the arc length from the conventional 
null charge surface, where the positive and negative values indicate outside and inside, 
respectively, with respect to the null surface. 

Combining FiglSland FiglHl we find a tendency that the field-aligned electric field becomes 
strong with increasing trans-field thickness. For example, in case45 2, case45 3 and case45 4, 
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in which jtot ~ 0.2, case45 3 has both the largest trans-field thickness and the strongest 
field-aligned electric field, on the other hand, case45 4 has both smallest trans-field thickness 
and the weakest field-aligned electric field. The available potential drop (= (5$a ) in the 
outer gap may be written as 5(^a ~ ^a^^gap/^p, where $a is the available potential drop 
on the stellar surface, 69gap is the angle of the foot points of the open field lines through 
the gap, and 6p is the angle of the polar cap. Because the angle SOgap increases with the 
trans-field thickness, the potential drop increases with the trans-field thickness. As shown 
in Paper2, the strength of the field-aligned electric field does not depend on the gap width 
along the field lines very much, because the effective charge density (j+ —j-)B/(3\^ + Bz and 
the transverse term — A_l$„co in the Poisson equation almost balance out in the gap f^ ^2.ip . 
As a results, the strength of the electric field tends to increase with the trans-field thickness. 
This effect is originally pointed out by the vacuum gap model of CHR86. 

From Figiniand FiglHl we also find that the field-aligned electric field of case45 3 is weaker 
than that of case45 1, although the trans-field thickness of case45 3 is larger. This refiects 
the screening effects of the particles on the field-aligned electric field. The dimensionless 
effective charge densities near the outer boundary for case45 1 and case45 3 are written as 
Peff = Bjgl [3\\ — Bz and pe// = B{ig + — B^, respectively. For case45 3, the external 

positrons assist in the screening near the outer boundary. 

As we have seen, the combination of the effects of the gap thickness and the screening 
of the particles determine the strength of the field-aligned electric field in the gap. 

3.1.4 'y-ray spectrum 

As we see in FigEl because most pairs are created near the inner boundary, the emission 
region for the inwardly 7-rays in the gap is restricted around the inner boundary, while the 
emission region for the outwardly 7-rays is whole outer gap. Furthermore, we expect that 
almost inwardly propagating 7-rays are converted into pairs by the magnetic or photon- 
photon pair creations when they pass close to the star. Therefore, we assume that the 
contribution of the inwardly 7-rays to the spectrum on the Earth is ignorable. 

FigEI shows the calculated curvature spectra of the outwardly propagating 7-rays for 
case45 l~case45 4. The solid lines in Fig represent the expected spectrum, which are com- 
posed of the gap emissions (dashed line) and the emissions (dotted line) from outside the 
gap. The crosses in Figd show the EGRET phase-averaged spectrum of the Vela pulsar. 
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Figure 7. Calculated 7-ray spectra of the curvature radiation for case45 l~case45 4 of the Vela pulsars. The solid line show 
the expected spectrum, which are composed of the inside (dashed line) and outside (dotted line) emissions. The crosses show 
the EGRET phase-averaged spectrum. The parameter A<^ adjusts the calculated flux to the observation. 



The adjusted azimuthal spread angle to the observed flux for each case is A0 ~ 25 radian 
(case45 1), ~ 12 radian (case45 2), ~ 6 radian (case45 3) and ~ 43 radian (case45 4). 

From Figl?! we find that the outside emissions (dotted lines) of case45 1 and case45 3 are 
fainter than the inside emissions (dashed lines). As the outer gap geometries in FigOlshow, 
the outer boundary of case45 1 or case45 3 is located near the light cylinder. Therefore, the 
outside emission region is smaller than the outer gap. For case45 2 and case45 4, on the other 
hand, the outside emissions contribute to the total spectrum below IGeV, because the size 
of the outside emission region is similar to or larger than the gap size. 

Combining FigEland FigQ we see that the spectral hardness reflects directly the strength 
of the field-aligned electric field. For example, the strongest (or weakest) electric field of 
case45 1 (or case45 4) in all cases produces the hardest (or softest) spectrum. We find that 
the spectrum becomes to be soft as the currents increase. Especially, by comparing the 
spectra of case45 1 and case45 4, we see that the spectrum sharply softens as the electron 
injection at the outer boundary increases. By comparing the spectra of case45 1 and case45 3, 
on the other hand, we also find that the increasing in the positronic component gradually 
softens the spectrum. 
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Figure 8. The solved outer gap (case25^ 1) of (ji, J2) = (0, 0) and jcrit ^ 0.2. The solid line in each panel shows the trans-field 
distribution of the internal current (upper left), the gap geometry (upper right), distribution of the field-aligned electric field 
(lower left) and the total spectrum (lower right). For comparison, we also display the gap geometry and the total spectrum of 
case45 1. 



3.2 Dependence on the internal current jg 

In the last section, we have investigated the dependence of the spectral hardness on the 
external components (ji, j2)- To clear the dependence of the spectrum on the critical internal 
component jcrit, we also calculate the case that ji = j2 = 10~^ and jcrit ~ 0.2 (denoted as 
case45^ 1) and compare the results with that of case45 1 (ji = j2 = 10^^, jcrit ~ 0.12). The 
results of case45^ 1 are summarized in Fig|Hl The upper left, upper right, lower left and lower 
right panels show the trans- field distribution of the internal component, the gap geometry, 
the distribution of the field-aligned electric field, and the total spectrum, respectively. For 
comparison, we also display the gap geometry and the total spectrum of case45 1 with dashed 
lines in FiglHl 

As FiglHl shows, the larger critical internal component, that is case^j^ 1, is predicted by 
a gap having the smaller trans-field thickness and longer width along the magnetic field 
lines. We expect that the gap of case^j^ 1 has a longer gap width W\\ than that of case45 1, 
because the internal component jg increases with the gap width. For case45 1, we adjust the 
positions of the outer and upper boundaries so as to reproduce a marginal gap structure 
having jcrit ~ 0.12. The field-aligned electric field changes it's sign in the gap due to the 
screening effect of the particles if we shift the outer boundary outward relative to that of 
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case45 1 without changing the upper boundary. As FiglHl shows, therefore, the gap of case^j^ 1 
has a smaller thickness and a longer width than those of case45 1. 

The decrease of the gap thickness causes a large transverse term A^$„co ~ —^nco/D\. 
As discussed in Paper2, the transverse term reduces the screening effects of the electrons, 
because it acts as a positive charge in the Poisson equation. The internal component jg can 
increase until the transverse term marginally sustain the negative charged super-GJ region, 
where |p| > \pgj\i in the gap. The choice of a large critical internal component tends to 
result in a smaller thickness. 

Because the smaller trans-field thickness gives the smaller field-aligned electric field and 
furthermore because the larger internal component jg gives the stronger screening effects for 
the field-aligned electric field, the resultant 7-ray spectrum of case^g^ 1 is softer than that of 
case45 1. 

We summarize the discussions in § ^3.11 and § ^3.21 The spectrum becomes soft as the 
current components, {jcrit, ji, j2), increase due to the screening effect on the field-aligned 
electric field. Especially, the increase in the electron component j2 at the outer boundary 
sharply softens the spectrum. On the other hand, the increase in the positronic component 
ji at the inner boundary gradually softens the spectrum. By increasing the positronic com- 
ponent ji, therefore, we can obtain a brighter emission with similar spectral hardness (see 
§§ ^H.H.2|1 . This dependence of the spectral feature on the current allow us to diagnose what 
current flows in the gap of the observed 7-ray pulsars. 



3.3 Application to the Vela pulsar 

In the present two-dimensional model, the gap width and thickness are linked so that we 
can discuss the gap structure in terms of the spectral characteristic; i.e. the cut-off energy 
and the flux. By comparing the model spectra with the observed spectrum, we can diagnose 
the free model parameters, that is, the current components (jcrit, Ji, J2) and the inclination 
angle Uinc- In this paper, we apply the theory to the Vela pulsar. 

In §§ §H.H.H we study that the EGRET spectrum of the Vela pulsar predicts the cases 
that have only a few electron injections (j2 ~ 0) at the outer boundary. In §§ §H.H.2[ we discuss 
parameter spaces of ainc and ji, which explain the EGRET spectrum, with jcru ~ 0.1, which 
is moderate value in the present electrodynamical model. 
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3. 3. 1 Flow pattern of the external current 

From Fig|7[ we find that the spectrum of case45 2, in which ji = j2 = 0.05, jcrit ~ 0.12 and 
«mc = 45°, becomes to be very soft as compared with the observation; the spectral cut-off 
energy is ~ 0.5GeV for case45 2 and ~ 3GeV for the observation. We also find that the 
adjusted spread angle for case45 2, A0 ~ 12 radian, is unrealistic. Starting from these facts, 
we can rule out any outer gap having a large particle injection at the outer boundary for 
the Vela pulsar, as follows. 

Firstly, we do not expect the cases having a larger current components than that of 
case45 2, because the spectrum softens as the current increases as discussed in § ^H.1I and 

Secondly, we do not expect the cases that the electronic component j2 is the same with 
that of case45 2, but the positronic component ji and/or the internal component jcrit are 
smaller than those of case45 2. This is because the adjusted azimuthal angle A0 increases 
with the decreasing ji and/or jcnt due to decrease of the outwardly migrating particles. 
Therefore, the unreality on the spread angle such as A0 = 12 radian of case45 2 grows with 
decreasing ji and/or jcnt- 

Finally, we can show that the hardness of the spectrum does not depend on the inclination 
angle very much. The increase in the inclination angle shifts the null point on the last open 
line toward the stellar surface. This shift increases the strength of the magnetic field in the 
outer gap, and also decreases the gap thickness with the mean-free path because the X-ray 
field becomes dense and because the collision of the inwardly propagating 7-rays approaches 
the head-on. The former effect tends to increase the strength of the field-aligned electric field, 
the latter to decrease it, the two effects compensate each other. As a results, the spectral 
hardness does not depend on the inclination angle very much. In fact, there is a tendency 
that the spectral hardness softens with the inclination angle for the large electron component 
at the outer boundary. Therefore, we can rule out any inclination angles with j2 ^ 0.05 for 
the Vela pulsar. 

As a result of the above three arguments, we can rule out any large electronic injection 
cases, i2 ^ 0.05, and we suggest that only a few electrons (j2 ~ 0.0) are injected into the 
gap at the outer boundary for the Vela pulsar. 

For j2 ~ 0, the decreasing of the gap thickness with increasing the inclination angle is 
more moderate than the case of the large external electron case, because the mean free path 
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Figure 9. The calculated outer gap geometry in wide parameter space 45°< ainc <75°. We use J2 = and jdrt ~ 0.1. The 
thin solid and dotted-dashed lines show the last-open line and the conventional null charge surface, respectively. 




Figure 10. The calculated total spectra with various inclination angles and positron components. The external electron and 
critical internal components arc je = and jcrit ~ 0.1. The calculated spectral hardness is explained with ji ^ 0.08 for 
ainc = 45°, ^ 0.10 for Uinc = 55°, ^ 0.15 for ai„c = 65° and ^ 0.20 for Oinc = 75°. 



of the outwardly 7-rays is very longer than the gap thickness. As a results, we can see that 
the strength of the field-aligned electric field and resultant hardness of the spectrum tend 
to increase with the inclination angle for j2 ~ 0. 
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3.3.2 Diagnosis the current and the inclination angle 

Following the discussion in ^^ ^3.3.11 we fix the external electron component j2 at zero for the 
Vela pulsar. For j2 ~ with the present model, the solutions disappear if the critical internal 
component jcrit exceed several ten percent of the Goldreich- Julian value. In this subsection, 
therefore, we adopt jcrit = 0.1 as a moderate value. With jcrit = 0.1 and j2 = 0, we diagnose 
other two model parameters, that is, the external positron ji and the inclination angle Uinc- 

FigEl and FigHUl shows the gap geometries and the spectra, respectively, for various 
positronic components and the inclination angles; upper left panel summarizes the results 
of ji = 0.01 (solid), 0.05 (dashed) and 0.08 (dashed) with aj„c = 45°, upper right panel 
ji = 0.01 (solid), 0.05 (dotted) and 0.1 (dashed) with ainc = 55°, lower left panel ji = 0.05 
(solid), 0.10 (dotted) and 0.15 (dashed) with a = 65° and lower right panel ji = 0.10 (solid), 
0.15 (dotted) and 0.20 (dashed) with ainc = 75°. By comparing the observed cut-off energy, 
we can restrict the positronic component ji such as ji ^ 0.08 for ainc = 45°, ^ 0.1 for 
55°, ^ 0.15 for 65°, and ^ 0.20 for 75°. The result of the marginal case for each inclination 
angle is shown with the dashed line. The values of the azimuthal spread angle A0 of the 
gap written in FiglHland FigUHlare adjusted to the observed flux. 

We notice that the increase in the inclination angle permits larger positronic component 
ji. For only a few electronic injections at the outer boundary j2 ~ 0, the field-aligned electric 
field in the gap tends to become strong with increasing the inclination angle as described 
in ^^ ^3.3. 1[ Stronger field-aligned electric field makes us possible to have larger positronic 
component ji, because we can increase ji until the particles screen the field-aligned electric 
field so that the calculated spectral cut-off energy does not explain the observation. 

For each positron external component ji of a^nc = 45° and ainc = 55° in FigEl the ad- 
justed azimuthal angles A0 are near or lager than 271 radian which are unrealistic. Therefore, 
we do not expect any cases of the critical internal component of jcrit ~ 0.1, which has been 
assumed in FiglHl for ainc = 45° and ainc = 55°. For smaller critical internal component 
than ji ~ 0.1, the unreality of the azimuthal angle increases. For more large critical internal 
component such as jcrit ~ 0.2 of case^^^ 1 in S ^3.2l the calculated spectra become to be soft 
as compared with the observation. On these ground, therefore, we safely rule out any cases 
of ainc ^ 55° for the Vela pulsar. 

For ainc = 65°, The adjusted azimuthal angle ~ 4 radian of the marginal case 
(dashed line) is slightly larger than A0 ~ vr. Because the three-dimensional gap geometry 
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affects the solid angle of the 7-ray beam and the expected flux on the Earth, we could deem 
the cases of ainc = 65° acceptable for the Vela pulsar in the present two-dimensional model. 
We expect that the large positron component such as ji ~ 0.1 runs through in the gap. 

For ainc — 75°, the adjusted azimuthal angle of Acf) ~ 2.5 radian for ji — 0.2 explains 
the observation very well. Furthermore, for ji — 0.1 (the sohd hue), the expected angle 
A(j) ~ 2.8 radian is also consistent with the observation. Therefore, we find that both the 
observed spectral cut-off energy and flux are explained by ji ^ 0.2 for ainc = 75°. As we 
have seen, the upper limit ji = 0.2 for ainc = 75° is determined by the spectral cut-off 
energy. On the other hand, the lower limit of the external positron ji is determined by a 
comparison between the model and the observed flux. Because the expected flux depends on 
the three-dimensional geometry, we should use a three-dimensional model to give the lower 
limit of ji for ainc = 75°. 

In summary, we have shown by solving the two-dimensional gap structure that the 
strength of the field-aligned electric field and the hardness of the produced 7-ray spectrum 
are strongly affected by the current through the gap. By comparing the model spectrum 
with the EGRET observation, we have successfully constrained the model parameters (the 
current and the inclination angle) for the Vela pulsar, as follows. Only a few external elec- 
trons (i.e. j2 ~ 0) injected into the gap at the outer boundary are predicted by the EGRET 
observation. We also safely rule out any cases of a^nc ^ 55°. With uncertainty due to the 
three dimensional effects, we conclude that the Vela pulsar has the inclination angle ainc 
greater than 65°. With the moderate critical internal component jcrit ~ 0.1, the positronic 
component ji is restricted such as ji ~ 0.1 for ainc — 65° and ji ^ 0.2 for ainc — 75°. 

We can also explain the observed spectrum with greater inclination angles. Following the 
RY95 arguments, the inclination angle will be not close to 90°, because the radio pulsation 
of the Vela pulsar shows a single peak. 

4 DISCUSSION 

4.1 Other high-energy emission processes 

In this section, we show the validity of the assumption that the synchrotron and the inverse- 
Compton process arc less important for the gap electrodynamics, and calculate the inverse- 
Compton spectrum to compare the result with the recent observation in TeV range. 
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4-1.1 Synchrotron process 

Pairs created by the pair creation process in the gap have significant pitch angle x ^i-nd will 
radiate the synchrotron photons. Initial energy of the pairs is several GeV, which is typically 
energy of the 7-ray photons radiated by the curvature process of the accelerated particles. 
After the birth, the pairs are accelerated to F ~ 10^'^ from the initial Lorentz factor of 
To ~ 10^'^ and will have pitch angle of x ~ To/r- It follows that the ratio between the 
synchrotron power and the curvature power becomes 

Psyn ^ f eBR,smx \' 3 ( ( J_Y' ( ( ^ ( J^V (a',) 

P,urv V rmec2 / ^ Uo^V Uo^V \WG ) UoOs-V WhRiJ ■ ^ ' 

This indicates that the synchrotron process is less important for the gap electrodynamics. 
However, it is not denied an another solution in which the synchrotron radiation play an 
important role. For example, as we have seen, for the large current case, the inner boundary 
is shifted toward stellar surface, where the synchrotron process would be important because 
the magnetic field strength significantly increases. The present model have predicted that 
the inner boundary of the gap that is favored for the Vela spectrum does not locate near 
the stellar surface (FiglHl). 

The large number of pairs will be produced outside the gap. These pars radiate syn- 
chrotron photons, the typical energy of which is 

= ^^M!^ sin X = 1.8 X 10^ ( ( sin x eV. (44) 

It is, therefore, entirely justified to neglect the synchrotron component emitted outside the 
gap when we considerer the 7-ray spectrum above 100 MeV. 



^.i.^ Inverse- Compton process 

TeV photons are emitted by the inverse Compton (IC) scatterings of infrared (IR) photons off 
the relativistic electrons and positrons, which emit GeV photons via the curvature process. 
For the Vela pulsar, Romani (1996) predicted the pulsed TeV flux at about 1% of the pulsed 
GeV flux, which has been ruled out by Air Cherenkov telescopes (Konopelko et al. 2005). 
Therefore, the theoretical prediction of TeV flux is important for checking the validity of 
the model. 

We assume the IR fleld is isotropic and homogeneous, because the IR photons will be 
radiated by the synchrotron process of the pairs with a larger pith angle. For isotropic 
case, the relativistic particle upscatters the soft photons to produce the following spectrum 
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(Blumenthal & Gould 1970) 

dN 3 c dNindEiR 

-aT- 



dtdE^ 4 r2 dEjR Ein 

2glng + (1 + 2g)(l - q) + %'^^}\~^^ ] , (45) 

where = ^En^/irieC^, q = £'i/re(l — Ei) and Ei = Ej/TemeC^ . From IR/Optical 
(Mignani & Caraveo 2001; Shibanov et al. 2003) and UV (Romani, Kargaltsev & Pavlov 
2005) observations, we adopt the following power low IR spectrum 

with Njfi = 1.3 X 10^^((i/0.25kpc)^. Integrating the upscattering over the high-energy par- 
ticles both inside and outside the gap, we obtain the TeV inverse Compton spectrum. 

Figim shows the expected inverse Compton spectrum with the curvature spectrum. The 
solid line and dotted line show the spectrum for (amc, ji) = (65°, 0.15) and (75°, 0.2), 
respectively, with {j2,jgap) = (0,0.1). We find that the model TeV flux less than 0.01% of 
the GeV flux is consistent with the recent upper limits by II.E.S.S.(Konopelko et al. 2005). In 
the present model, the high-energy particles accelerated in the gap make both the curvature 
and inverse- Compton spectra. For such a case, the flux ratio between the inverse-Compton 
and curvature processes does not depend on the inclination angle and the gap size, because 
the ratio is determined by ratio between emission rates of a particle so that 
Lie ^ rm^c'aTEin^ 

T ~ 2 e^r-t 

J^curv 

= 6.4 X 10- lE,„<lN,,liE,„\ 

\W) \wiJ \ lOiicm-3 j ^ ^ 

It is clear that the inverse-Compton process is less important for the gap electrodynamics. 



4.2 Comparison with the previous works 

We have estimated the inclination angle that aj„c ^ 65° for the Vela pulsar. This result 
is marginally consistent with the result by RY95, in which ainc ~ 65° for the Vela pulsar 
was predicted to explain the phase difference between the radio and 7-ray pulses. RY95 
considered the gap having no external component (ji = j2 = 0) and starting from the 
conventional null surface. However, the present model has predicted the particle injection 
at the inner boundary (ji 7^ 0,j2 = 0) to explain the observed flux, although the obtained 
inner boundary is located near the conventional null surface (FigEI). 
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Figure 11. Pulsed gamma-ray spectrum from the Vela pulsar. The solid line and dotted line show the spectrum for (ai,ic, jl) = 
(65°, 0.15) and (75°, 0.2), respectively, with (j2,igap) = (0,0.1). The upper limits in 0.1-lTeV range has been determined by 
H.E.S.S. (Konopelko et al. 2005). 

RY95 have assumed that the gap is extending to the hght cyhnder to produce the first 
peak of the light curves. In this paper, the gap geometry in FiglHlthat is favored for matching 
the spectrum of the Vela pulsar does not predict such long outer gap. However, we note that 
the accelerated particles in the outer gap emit the 7-ray photons outside the gap. As Fig|7| 
shows, the curvature photons emitted outside the gap by the outwardly moving particles also 
contribute to emissions above lOOMeV. Specifically, the outwardly moving particles escape 
from the light cylinder with the Lorentz factor F ~ 1.2 x 10'^(i7/100s~^)~^/^(< Rc > jwi^'^l'^, 
the corresponding curvature photon energy is E^^ ~ 160(< Rc > /voi^ MeV, where < Rc > 
is the averaged curvature radius of the field line, along which the particles move. Therefore, 
although the outer boundary of the gap is located inside of the light cylinder, the high-energy 
photons are emitted near the light cylinder and the double peak structure in the calculated 
light curves is expected. The detail calculation will be done in the subsequent papers. 

Dyks & Rudak (2003) proposed the radiation emission region starting from near the stel- 
lar surface to explain the presence of the outer-wing emissions and of the off-pulse emissions 
of the Crab pulsar. With the two-dimensional outer gap model, we found that such gap has 
a large particle injection at the outer boundary (Paper2). For such cases, ji ^ j2, however, 
the our model (e.g. case45 4) also predicts the softer spectrum than the observed one of the 
Vela pulsar, as we discussed in ^ ^3.31 
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4.3 Screening of the electric field 

As we have seen, the inner boundary, where E\\ = 0, is located at the position where the 
local charge density caused by the current is equal to the Goldreich- Julian charge density, 
that is, (ji — j2 + jg)B ~ Bz is satisfied on the inner boundary. However, the condition 
I ^ (ji ~ J2 + jg)B\ < \ — Bz\ is satisfied between the stellar surface and the inner boundary. 
This negative charge depletion from the Goldreich- Julian charge density predicts a field- 
aligned electric field between the stellar surface, on which $„co = is defined, and the 
inner boundary. In such a case, there is a potential drop along the magnetic field lines. In 
the present paper, however, we have imposed the condition (jSZ|), ^nco = 0, at the inner 
boundary with the assumption that a screening mechanism of the electric field is working 
and that the potential drop between the stellar surface and the outer gap is much smaller 
than the potential drop in whole outer gap. In this subsection, therefore, we discuss the 
screening mechanism between the stellar surface and the inner boundary with the case that 
ji = J2 = 0. 

Firstly, as the particles that screen the electric field, we mention a population of the 
electrons trapped between the stellar surface and the outer gap, such as the electron clouds 
trapped between the stellar surface and the boundary of the outer magnetospheric gap 
appeared in Smith, Michel & Thacker (2001). These electrons stabilize the force free surface 
(i.e. the inner boundary) between the plasma side and the outer gap side. Therefore, we 
expect that three kind of the particles will populate between the star and the inner boundary; 
the primary particles accelerated in the gap, the pair particles produced outside gap by the 
inwardly 7-rays and the trapped electrons. The vacuum outer gap model of CHR86 also 
assumed these trapped electrons between the stellar surface and the traditional null charge 
surface. Such trapped electrons will be distributed to fill the negative charge depletion, 
jgB — Bz, from the Goldreich- Julian charge density, and to satisfy the condition = 
between the stellar surface and the inner boundary. 

Because the inner boundary will be affected by the plasma condition around the gap, we 
consider a particular case that there are no trapped electrons as an another possibility, that 
is, only primary and pair particles are migrating toward the stellar surface along the magnetic 
field line linking the star and the outer gap. With the assumption that the screening distance 
along the field line is much smaller than the trans-field thickness, we take the field-aligned 
part in the Poisson equation to obtain the field-aligned electric field. In this case, we can 
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produce the solution that the field-ahgned electric field E\\ undergoes spatial oscillation, and 

therefore is screened out by the small velocity difference of the pair particles. Because the 

amplitude of the potential (~ 10^ in the non-dimensional value) is much smaller than the 

potential drop (~ 10^) in the outer gap, the condition (|H7j). ^nco = 0, at the inner boundary 

is a good treatment. Following Shibata, Miyazaki & Takahara (2002), however, we can show 

that a surprisingly large pair-creation rate is required to screen. In fact, the screening is 

achieved if the pairs were created with a rate 
/ dm \ ( f2 \ / \ 2 



> 0.0025 f " .) /7<,y/al^^ I 

VlOOrads^V VlOV \ r ) ' ' 



where m is the pair flux divided by ^Binj'^'^ ■, -Bm is the strength of the magnetic fleld at 
the inner boundary and 70 is the initial Lorentz factor of the pair. This rate indicates that 
the total number of the pairs produced within the typical screening distance ~ lO^cm 
is Np/Ncj ~ 10^, where Nqj represents the typical Goldreich- Julian value. Such a very 
large pair creation rate has not been predicted by the present cascade model, in which 
dm/ds\\ ~ lO"*^ /cm (FigEJ. On these ground, we expect that the pair polarization does not 
screen out the fleld-aligned electric fleld within the small distance along the fleld line. 

As discussed above, if no trapped electrons exist between the stellar surface and the 
outer gap, the fleld-aligned electric fleld will be caused outside the gap over a wide range 
along the fleld line. For such a case, the trans-fleld part in the Poisson equation determines 
the non-corotational potential, and the fleld-aligned part is not important, so that 

[-PO+P+ + P-], (48) 

where po = jgB — B^, p+ and p_ are, respectively, the charge densities of the par positrons 
and electrons, and is the distance in the trans-fleld direction. By ignoring the trans-fleld 
distribution of the charge density, the solution of equation becomes $nco ~ [— Po + P+ + 
P-]C,±iD± — ^±)/2, where we use the boundary conditions that $„co = at = and 
^± = D^, where D± represents the trans-fleld thickness. The fleld-aligned electric fleld is 

^\\ = ^- ^(-Po + P+ + P-) ^ , (49) 

where .^n is the fleld-aligned distance increasing toward the star. 

Let us consider the sign of the fleld-aligned electric fleld of equation ()49|1 . If we ignore the 
contribution of the pair particles on the electric fleld, we flnd that the electric fleld becomes 
a positive because the outer gap is located on the fleld lines that curve away from the 
rotational axis and because the derivative —dpo/dC,\\ has a positive value between the stellar 
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surface and the inner boundary. This positive electric field accelerates the pair electrons and 
decelerates the pair positrons, which are migrating toward the star. The continuity equations 
of the pair particles predict the field-aligned distribution of the charge density increasing 
toward the stellar surface, that is c?(p+ + pJ)/d^\\ > 0. The derivative of the total charge 
density in the right hand side of equation (jl^ has positive value, and therefore the effect of 
the pair polarization strengthens the positive electric field, which will appear unless there 
are such pair particles. Although some positrons are returned and are accelerated toward the 
gap by the positive electric field, the positive electric field will be held between the stellar 
surface and the inner boundary even if the effect of the discharge of the pairs is taken into 
account, as has been indicated by Hirotani (2005, in preparation). 

Recently, Hirotani (2005) has solved the structure of the particular outer gap with the 
particle motion by setting the inner boundary at the stellar surface. He has indicated that 
such outer gap is composed of two region; one region (called "minor part") takes on a 
minor part of the potential drop of whole outer gap and has a small positive field-aligned 
electric field, and other (called "major part") takes on a major part of the potential drop 
and has a large field-aligned electric field. The major part, which corresponds to so-called 
outer gap in the present paper, is located in the outer magnetosphere and the minor part 
is extending between the stellar surface and the major part. Because the potential drop in 
the minor part of the gap is negligibly small in comparison with that in the major part, the 
boundary condition (j37j) . $„co = pTj) . at the inner boundary of the major gap will be a 
good treatment as long as we are concerned with the major part of the gap. 

With above discussions, therefore, we conclude that the boundary condition that $„co = 
at the inner boundary, where = 0, for the outer gap is a good treatment. 



4.4 The connection between the global and local models 

As we have seen, the outer gap structure and produced 7-ray spectrum are controlled by the 
current, which circulates in the magnetosphere globally. By comparing the observation, our 
local model can constrain the current fiow pattern in the gap and the inclination angle, which 
cannot be observed directly. For further discussions, we would need a connection between 
the global and local models, as follows. 

In the previous sections, we had considered the outer boundary following a curve line 
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Figure 12. Outer gap structure with (j'l, J2)=(0.1,0) and ainc = 65°. Each panel is same with Fig|S] The sohd lines show 
the results of the gap having the outer boundary located near the light cylinder at the lower part. The dashed lines show the 
results of the gap having the outer boundary of (^out = 0.64. 



perpendicular to the magnetic field lines, that is C,out =constant. It is notable that the shape 
of the outer boundary is controlled by the trans-field distribution of the total current. 

Around the lower boundary, only a few pairs are produced by the pair creation, because 
the 7-rays are always emitted to the convex side of the magnetic field lines. Therefore, the 
shape of the outer boundary can become such as the gap structure displayed with solid lines 
in FiglT^ of which the outer boundary changes from C^out = 0.64 to (out = 1 in region within 
16% of gap thickness measured from the lower boundary. For comparison, we show the gap 
structure (dotted line) having the outer boundary of (out = 0.64. 

The electrons produced in the extending region emit the inwardly propagating 7-rays to 
convex side of the magnetic field lines and make pairs around the inner boundary. In the 
upper left panel of Fig lT^ we can see that an effect of these new born pairs is the increase 
of the internal current around the hight {x — Xi)/{Xup — Xi) ~ 0.2 measured from the lower 
boundary. FiglT^ The difference between the two currents produces the difference between 
the outer boundaries in the upper right panel. 

We find in Fig^J that the calculated spectra for both cases are consistent with the 
observation. Therefore, it is difficult to deduce the shape of the outer boundary using the 
present local model. As a global study, CHR86 postulated that the outer gap is extending 
until the light cylinder, such as the gap geometry depicted with the solid line in Fig lT^ 
due to the presence of an ionic (or a positronic) particle outflow from the light cylinder. 
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For the aligncd-rotator, on the other hand, Mestel et aL (1987) studied a structure of the 
magnetosphere with a global current and an acceleration region, which closes within the 
Ught cyUnder (see also Mestel 1999). These global studies solve the global current, which 
interacts with the polar cap, the outer gap and the pulsar wind. 

Because the local model can predict the current through the gap with the observations 
as we have seen in the present paper, and because the global study can solve the global 
current flow in the magnetosphere, the combination between the local and global studies 
would give a consistent understanding of the structure of the active pulsar magnetosphere. 
For such a discussion, a three-dimensional electrodynamical model is needed to compare 
with the observed phase-resolved spectra and light curves. 
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